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ABSTRACT 

X-ray flares were discovered in the afterglow phase of gamma-ray bursts (GRBs) 
by the Swift satellite a decade ago and known as a canonical component in GRB X-ray 
afterglows. In this paper, we constrain the Lorentz factors of GRB X-ray flares using 
two different methods. For the first method, we estimate the lower limit on the bulk 
Lorentz factor with the flare duration and jet break time. In the second method, the 
upper limit on the Lorentz factor is derived by assuming that the X-ray flare jet has 
undergone saturated acceleration. We also re-estimate the initial Lorentz factor with 
GRB afterglow onsets, and find the coefficient of the theoretical Lorentz factor is 1.67 
rather than the commonly used 2 for interstellar medium (ISM) and 1.44 for the wind 
case. We find that the correlation between the limited Lorentz factor and the isotropic 
radiation energy of X-ray flares in the ISM case is more consistent with that of prompt 
emission than the wind case in a statistical sense. For a comparison, the lower limit on 
Lorentz factor is statistically larger than the extrapolation from prompt bursts in the 
wind case. Our results indicate that X-ray flares and prompt bursts are produced by 
the same physical mechanism. 

Subject headings: gamma ray: bursts - radiation mechanism: non-thermal 


1. Introduction 

Gamma-ray bursts (GRBs) are the most luminous explosive events in the Universe. It is 
well known that the observed gamma-ray emission are produced by relativistic outflows. The 
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fireball of a GRB is required to have a relativistic speed toward the earth in order to avoid the 
“compactness problem”. The Lorentz factor of the fireball increases with radius in the radiation- 
dominated acceleration phase. At the end of this phase, the fireball enters the matter-dominated 
coasting phase, which maintains a constant Lorentz factor (for reviews see Piran 1999; Zhang 
& Meszaros 2004; Meszaros 2006; Kumar & Zhang 2015). The Lorentz factor in the matter- 
dominated coasting phase, called the initial Lorentz factor, is a crucial parameter in understanding 
the physics of GRBs. There are several methods to constrain this initial Lorentz factor T 0 , but 
the most effective method is that the peak time of a GRB early afterglow onset is taken as the 
deceleration time of the outflow. An estimation of this initial value is possible by measuring the 
peak time of the afterglow light curve and the prompt isotropic energy E jiso (e.g., Molinari et al. 
2007). This method was successfully applied in Liang et al. (2010) with a large sample of the early 
afterglow light curves. Liang et al. (2010) discovered a tight correlation between T 0 and E yiso . In 
this paper, we re-estimate the initial Lorentz factor of GRBs in the ISM and wind cases, and also 
obtain a tight correlation between T 0 and E l iso . The T 0 and E ltiso correlation in the wind case is 
even more tighter than that of the ISM case. 

X-ray flares are common phenomena in GRB X-ray afterglows in the Swift era. According 
to the temporal behavior and spectral properties of X-ray flares, it is widely argued that X-ray 
flares are produced by long-lasting central engine activities (Burrows et al. 2005; Ioka et al. 2005; 
Fan & Wei 2005; Falcone et al. 2006; Zhang et al. 2006; Wang & Dai 2013). It is generally 
believed that X-ray flares have the same physical origin as the prompt emission of GRBs. X-ray 
flares usually happen at ~ 10 2 - 10 5 s after prompt emissions, transferring a large percentage of 
the outflow energy to the radiation. Most of the methods proposed to constrain the initial Lorentz 
factor of GRBs are not suitable for X-ray flares. Since the fluence of most X-ray flares are smaller 
than the prompt emission, their energies and Lorentz factors are supposed to be smaller than those 
of GRBs. In this paper, we use two different methods to place limits on the Lorenz factor of X-ray 
flares. In Section 2, we introduce the method to estimate of the initial Lorentz factor r 0 of GRBs. 
The two methods of constraining the X-ray flares Lorentz factors are shown in Section 3. Sample 
study are presented in Section 4, and our results are summarized and discussed in Section 5. A 
concordance cosmology with parameters Ho = 71 km s _1 Mpc -1 , ELm = 0.30, and = 0.70 are 
adopted. Notation Q n denotes <2/10" in the cgs units throughout the paper. 


2. Estimating the Initial Lorentz Factor r 0 

When a relativistic fireball shell sweeps up the circumburst medium, two shocks appear: a 
reverse shock propagating into the fireball shell, and a forward shock propagating into the ambient 
medium. We suppose that the fireball shell and two shocks are spherical. Physical quantities 
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denoted by “/” are defined in the co-moving frame. We divide this two-shock system into 4 regions 
(Sari & Piran 1995; Yi et al. 2013): (1) The unshocked ambient medium (n\, e\, pi, TO, (2) the 
shocked ambient medium (. n 2 , e 2 , p' 2 , T 2 ), (3) the shocked fireball shell (n 3 , e' 3 , p 3 , T 3 ), and (4) 
the unshocked fireball shell (n A , T 4 = T 0 ), where n is the number density, e is the internal energy 
density, p is the pressure, and T is the bulk Lorentz factor. Hereafter we calculate the properties of 
the reverse shock emission at the radius R x , where the reverse shock finishes crossing the ejecta 
shell and the Lorentz factor of the shell is T x . 

The assumptions of the equilibrium of pressures and equality of velocities along the contact 
discontinuity lead to p 2 = p 3 and T 2 = T 3 , respectively. With the jump condition for the shocks and 
the equilibrium of pressures, we can obtain, 


(4T34 + 3XT34- 1)«4 - (4T 2 + 3)(r 2 - l)«i. 


The Lorentz factor of the reverse shock T 34 could be approximated as 


r 34 = 


1 

2 



i/r|+rj\ 

2 v r 4 r 3 ) ’ 


(1) 

( 2 ) 


as long as T 3 S> 1 and T 4 S> 1. Substituting Equation (2) into Equation (1), we can obtain the 
following equation 


(r 4 -r 3 )- 

r 4 r 3 


cr 4 +r 3 r 

r 4 r 3 


n A = ATlni. 


( 3 ) 


Because T 3 1 , T 4 S> 1 and T 4 > T 3 , we ignore the constant 1/2 term in the Equation (3) and 
thus we can obtain the solution of this equation (ignoring the negative solution, also see Panaitescu 
& Kumar 2004) 

r 4 


1+2 r 4 (fti/n 4 ) 


1 / 2 ' 


1 / 2 ' 


( 4 ) 


Here we obtain the relation between the Lorentz factor of shocked fireball shell T 3 and the initial 
Lorentz factor T 0 (T 4 ), which depends on the ratio of these two comoving densities. The number 
density of the ambient medium is assumed to be ri\ =AR~ k (Dai & Lu 1998; Meszaros et al 1998; 
Chevalier & Li 2000; Wu et al 2003, 2005; Yi et al. 2013), such a circumburst medium is a 
homogeneous interstellar medium (ISM) for k = 0, and a typical stellar wind environment for k = 2. 
The fireball shell is characterized by an initial kinetic energy £)., initial Lorentz factor T 4 , and a 
width A in the lab frame attached to the explosion center, so the number density of the shell in 
the comoving frame is n A = E k /(47im p c 2 R 2 ATf). The ratio of the comoving number density of the 
relativistic shell n' 4 to the number density of the ambient medium n A defined in Sari & Piran (1995) 
is 

__ X 

n j AnAnipC 2 AT \R 2 ~ k A T 2 A R 2 ~ k ’ 


( 5 ) 



- 4 - 


where X = Ek/( 47 rAm p c 2 ). The difference between the lab frame speed of the unshocked fireball 
shell and that of the reverse shock is (Kumar & Panaitescu 2003), 


lAfTlmV 1 - 4/13 5 


Pa-Prs-^t —— -^7 ' 


I 4 V < 


r 4 \f 


( 6 ) 


Considering the thin shell case A ~ R/{2T\), we can calculate the radius R x where the reverse 
shock finishes crossing the fireball shell. 


A (R x ) = 


(P 4 - /3 RS )dR. 


(7) 


/ 0 


The substitution of Equation (5) and (6) into Equation (7) leads to 

[2(3-k) z Xl ~ 
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So the comoving density ratio at R x is 


/x=^ = 


2X 


5.6- 
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(9) 


ni R 3 x k (5 -k ) 2 

Substituting Equation (9) into Equation (4), we can obtain the Lorentz factor of the reverse shock 
as it finishes crossing the shell 

r 4 r 4 


r x = 


1/2 


[1+0.357(5-k)] 


1+2T 4 (/C-7(2X)) 1/2 
Therefore, the relation between T x and the initial Lorentz factor T 0 is 

T x =0.60T 0 , for k = 0(ISM), 

and 
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T x = 0.70To, for k = 2 (Wind). 

For the thin shell case, the reverse shock crossing time T x is almost corresponding to the deceler¬ 
ation time Tdec, .i.e, T x ~ T dec . Therefore, we can derive the initial Lorentz factor in the ISM and 
wind type cases (also see Panaitescu & Kumar 2004). For k = 0 (ISM), 


T 0 = 


1 


0.60 


1.67 


3 E ■ 


and for k = 2 (wind), 


Tn = 


0.70 


:T X = 1.44 
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With the isotropic-equivalent energy E lpso and the peak time of the afterglow onset t PjZ , we can 
estimate the initial Lorentz factor of GRBs, where t pz = t p /{\+z). Liang et al. (2010) discovered a 
tight correlation between r 0 and E piso using 20 GRBs which show deceleration feature in the early 
afterglow light curves. Other work also confirmed this correlation, but with different methods and 
power-law indices (Ghirlanda et al. 2012; Lii et al. 2012). Using the data of t p - and E~ piso from 
Liang et al. (2010, 2013) and Lii et al. (2012), we re-constrain the initial Lorentz factor, and also 
discover a tight r 0 and E yiso correlation for the ISM and wind cases. The T 0 and E piso correlation 
in the wind case is even more tighter than that in the ISM case, as shown in Figs. 4 and 5. 


3. Methods of Constraining the Bulk Lorentz Factor of X-ray flares 


Because most of the models estimating the bulk Lorentz factor during the GRB prompt emis¬ 
sion phase are inapplicable for the X-ray flares, we introduce two methods to constrain the upper 
and lower limits on the Lorentz factor of X-ray flares in this section. In principle, GRB outflows 
could be structured, but we suppose that the outflows are conical (also called jet-like) and the 
half-opening angles of the outflows have a same/constant value in a single GRB, i.e., the half¬ 
opening angle of the X-ray flare jet is the same as that of the prompt jet in one GRB. Therefore, 
we suppose that each X-ray flare is produced in a conical uniform jet with half-opening angle 9 j. 
If there are several flares in one X-ray afterglow, we also assume these X-ray flares have the same 
jet half-opening angle 0j. 


Method I: Lower limit on the Lorentz factor of X-ray flares. We use the late internal shock 
emission model to constrain the lower limits of the Lorentz factor of X-ray flares (Wu et al. 2006, 
2007). The quick decline of X-ray flares after the peak time is widely interpreted as the high 
latitude component of X-ray pulses (Burrows et al. 2005; Tagliaferri et al. 2005; Zhang et al. 
2006; Nousek et al. 2006; Liang et al. 2006). Wu et al. (2007) supposed that emission from 
the same internal shock radius Rj nt but with different angles would have different arrival times to 
the observer, which is due to light propagation effect. The delay time of different photons with 
different 9 {6 < 9j) is 


AT - 


R int (l - COS#) 


(15) 


c 

If photons are emitted at an angle 9 = 9j , then the delay time AT is about the duration time of the 
X-ray flare. We denote the timescale of the decay part of an X-ray flare by T decay . We could put 
forth a constraint on the decaying timescale of the flare and the jet-opening angle, i.e. 


Tdecay *7 


R int (l-cos6>) 


c 


(16) 
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Another timescale is the angular spreading variability time T rise , which is the rise time of an X-ray 
flare, 


T = 

± rise 


R int (l-cos(l/T*)) R m 


2ry 


(17) 


Combining the decaying timescale and rising timescale, we could get a lower limit on the Lorentz 
factor of each X-ray flare, 


r v > 


Tdecay \ 
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With the jet opening angle 9j, decaying timescale and rising timescale, we can constrain the lower 
limit of the Lorentz factor of each X-ray flare. The jet half-opening angle 6j can be estimated from 
the late afterglow of a GRB if there is a jet break in the afterglow light curve. So, we selected GRBs 
which have jet breaks and flares in the X-ray afterglow light curve in our sample. The sample is 
listed in Table 1. 


Method II: Upper limit on the Lorentz factor of X-ray flares. The physical mechanism of X- 
ray flares is still unclear. X-ray flares are known to be similar to those of prompt emission pulses 
through studying of the temporal behavior and energy spectrum (Chincarini et al. 2010). We gen¬ 
erally consider X-ray flares having the same physical origin as the prompt emission of GRB, and 
they are all due to the long-lasting activity of the central engine (Fan & Wei 2005; Zhang et al. 
2006). But there is still controversy about the origin of X-ray flares of GRBs. X-ray flares of short 
GRBs can be produced by differentially rotating, millisecond pulsars from the mergers of binary 
neutron stars (Dai et al. 2006). Magnetic reconnection-driven explosions lead to multiple X-ray 
flares minutes after the prompt GRB. Wang & Dai (2013) performed a statistical study of X-ray 
flares for long and short GRBs, and found energy, duration, and wait-time distributions similar to 
those of solar flares, which indicates that X-ray flares of GRBs may be powered by magnetic re¬ 
connection. According to the standard GRB fireball model, after the initial radiation-dominated ac¬ 
celeration phase, the fireball enters the matter-dominated “coasting” phase (Piran 1999; Meszaros 
2002; Zhang & Meszaros 2004). Whether the fireball is baryon-rich or not determines how long is 
the initial radiation-dominated acceleration phase. If the fireball is baryon-poor, the initial energy 
of the fireball will be quickly converted into radiation energy and produce bright and brief thermal 
emission, which is inconsistent with most of the observations. The spectra of X-ray flares are typ¬ 
ically non-thermal, with a photon index of about ~ -2.0 (Falcone et al. 2007). This suggests that 
X-ray flares happen when the jet is optically thin. Meanwhile, thermal emission just before any 
X-ray flare with compatible flux has never been detected, indicating that the jet responsible for the 
flare attains saturated acceleration. This requires baryon loading in the jet to be large enough, or 
that the Lorentz factor has an upper limit. On the other hand, the upper limit on the Lorentz factor 
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of an X-ray flare can be estimated as (Jin et al. 2010) 


F, < 


Lo 7 


Snm p c 3 Ro 


(19) 


which depends on the total luminosity L and initial radius Rq of the flare outflow. Jin et al. (2010) 
assumed that the observed X-ray flare luminosity is just a fraction (e x = 0.1) of the total luminosity 
of the outflow, that is L x = 0.1L (also see. Fan & Piran 2006). R 0 is taken to be 10 7 cm, which is 
comparable to the radius of a neutron star. R 0 = 10 7 cm is a conservative value, in some cases, R 0 is 
taken to be ~ 10 8 cm or even larger (Pe’er et al. 2007). With the proper values, we are able to get 
an upper limit on T x . In our sample, L x is taken as E X IW /T )() X , where E x iso and T. )(} x are the isotropic 
energy and duration time of one flare respectively. The redshift of those GRBs in our sample are 
all measured, the isotropic 0.3-10 keV energy of the X-ray flares in the sample can be estimated 
from the fluence as 


F = 


1 +z 


( 20 ) 


where S x is the fluence of an X-ray flare. 


4. Case Studies 

Our method for placing lower limits on the Lorentz factor is feasible if the X-ray light curve 
presents flare and jet break simultaneously. Our sample consists of 20 GRBs with X-ray flares, 
redshift and jet break time (Falcone et al. 2007; Chincarini et al. 2010; Bernardini et al. 2011; 
Lu et al. 2012). Some of them have several flares in one GRB. The total number of X-ray flares 
is 43. We assume that the opening angle is the same for the jets responsible for prompt emission 
and late X-ray flares in a single GRB. According to the appearance time of X-ray flares, we mark 
the corresponding numerical order, which can be seen in the Table 1. The T rise , T decay , S x , and T )i) x 
of each flare are reported in Falcone et al. (2007), Chincarini et al. (2010) and Bernardini et al. 
(2011). 9j md and 0™ are calculated with the data taken from Lu et al. (2012) when the Eq. (22) 
and Eq. (23) are applied. The Lorentz factor in the wind type circumburst media is (Chevalier & 
Li 2000) 

r “ 5,9 ( i r) ' <21) 

where E k 52 is the initial kinetic energy of the fireball shell in units of 10 52 ergs, t days is the observer’s 
time in units of days, A = M w /4ttV w = 5 x 10 n A* g cm -1 is the wind parameter, M w is the wind 
mass-loss rate, and V w is the wind velocity. Because the jet break effects are considered when 



-8 


T ~ 0/, so Equation (21) could be used to estimate the jet half-opening angle for the wind case, 

8*'"“' =0.12rad <“> 

where rj is the efficiency of prompt GRBs and A* = 1 adopted in this paper. The jet half-opening 
angle in the interstellar medium case could be described by (Sari et al. 1999; Rhoads 1999; Frail 
et al. 2001; Yiet al. 2015), 


ef M = 0.076 rad 


-2L) 

1 day / 


3/8 


l+z\ ^ 77 - 1/8 ( V \ V8 ( U \ 1/8 

2 J 'y, iso, 53 y 0 2 / Ucm-V 


(23) 


The distribution of jet half-opening angles for the ISM and wind cases is shown in Fig 1. 

We assume that X-ray flares are coming from relativistic jets. The tail of an X-ray flare is 
interpreted as emission from high latitude areas of the jet. The duration of the flare is determined 
by the half-opening angle of the jet through the curvature effect. With the jet half-opening angle 
6j estimated from the jet break time, decaying timescale and rising timescale of the flare, we can 
obtain the lower limit on the Forentz factor of the flare via Equation (18). The upper limit on the 
Forentz factor is determined by the total luminosity and initial radius of the outflow. The observed 
average luminosity is obtained from the isotropic 0.3 - 10 keV energy of the X-ray flare averaged 
by the duration time of the flare. The total luminosity of the flare is assumed to be 10 times of the 
observed X-ray luminosity. As already mentioned, the initial radius of the outflow R 0 is taken as 
10 7 cm. The obtained limits on the bulk Forentz factor of X-ray flares range from tens to hundreds, 
as can be seen in Fig 2 and 3. We find that in the ISM case the correlation between the Forentz 
factor and the isotropic radiation energy of X-ray flares is almost consistent with that of prompt 
emission of GRBs (Fig 4). However, in the wind case the lower limit on the Forentz factor is 
statistically larger than the extrapolation from prompt bursts (Fig 5). 


5. Discussion 

X-ray flares are common features in GRB X-ray afterglows, and most of them have occurred 
at the early period. We can conclude that all the flares in our sample occurred before the jet break, 
which can be seen from the Table 1. Here, we define / = 0j hl /0j hx , where 6j iX is the half-opening 
angle of the X-ray flare jet while Q j hl is the half-opening angle of the jet responsible for prompt 
emission. If / = 1, one suggests that the jet may be conical and the flare jet and prompt emission 
jet have the same half-opening angle as we discussed above. The jet opening angle might be larger 
during the prompt emission and smaller for the X-ray flares, i.e., / > 1, which has been predicted in 
some models with magnetic-dominated jets (Fevinson & Begelman 2013; Bromberg et al. 2014). 
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In this case, the lower limit of flare Lorentz factor would be larger than that estimated with Eq. 
(18) assuming d j x = 8j n . The corresponding lower limits on X-ray flare Lorentz factor in Figs. 4 
and 5 will increase by a factor of /, making the X-ray flares possibly more inconsistent (especially 
for the wind case) with the extrapolation of the correlation between isotropic radiation energy and 
Lorentz factor of prompt emission of GRBs. On the other hand, although several observations 
suggest that in some GRBs the ejecta may have large scale magnetic fields and therefore the ejecta 
could be magnetized, the degree of magnetization is usually estimated as a < a few in the afterglow 
phase. So in this paper we assume that for simplicity GRB jets have negligible magnetization, and 
the outflows have same half-opening angle (/ = 1) in one GRB. 

In addition, the fluence of most X-ray flares are smaller than that of prompt emission, their 
energies and Lorentz factors are supposed to be smaller than those of GRBs. The initial Lorentz 
factor of GRBs in this paper is generally larger than a few hundreds, and it is always larger than 
the lower limits of X-ray flare Lorentz factor in the same GRB. In Fig 6, we plot 5 GRBs having 
prompt and flare Lorentz factors, GRBs 050820A, 060418, 060906, 070318, and 071031. The 
initial Lorentz factors of these 5 GRBs are generally much larger than the lower limits of the 
Lorentz factor of their X-ray flares, and usually smaller than the upper limits of flare Lorentz 
factors. 


6. Conclusion 

The initial Lorentz factor is a key parameter to understanding the GRB physics. In this paper, 
we have re-estimated the initial Lorentz factor in a more accurate way. From Equation (13), we 
obtain a coefficient 1.67 for the ISM case, instead of 2 adopted in previous literature. We also 
constrain the initial Lorentz factor in the wind case, which is shown as Equation (14). With the 
estimated initial Lorentz factors in this paper, we confirm the tight correlation between the initial 
Lorentz factor and isotropic energy of GRBs for the ISM case. There is an even tighter correlation 
between the initial Lorentz factor and isotropic energy of GRBs for the wind case. 

Our sample consists of 20 GRBs with X-ray flares, whose redshifts and jet break times are 
known. Some of them have several flares. The total number of X-ray flares in our sample is 43. 
We assume that the half-opening angle is the same for the jets responsible for prompt emission and 
late X-ray flares in one GRB. Our results are shown in Fig 4 (ISM) and Fig 5 (Wind), which also 
show the correlation between isotropic radiation energy and the Lorentz factor of prompt emission 
of GRBs. The obtained limits on the bulk Lorentz factor of X-ray flares range from a few tens to 
hundreds, together with the isotropic radiation energy, are generally consistent with the correlation 
for prompt GRBs for the ISM case. Our results indicate that X-ray flares and prompt bursts may 
be caused by the same mechanism, as both are produced by the long-lasting activity of the central 



engine. However, in the wind case the lower limit on Lorentz factor is statistically larger than the 
extrapolation from prompt bursts. 
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Fig. 1.— The distribution of half-opening angles of GRB jets. The solid and dash lines are corre¬ 
sponding to the ISM and wind cases, respectively. 



log (lower r) 


Fig. 2.— The distribution of the lower limits on the Lorentz factor of X-ray flares for two types of 
circumburst environment. The lower limits on Y x range from tens to a few hundreds for the ISM 
(solid) and wind (dashed) cases. 
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Fig. 3.— The distribution of the upper limits on the X-ray flare Lorentz factor. 



Fig. 4.— Lorentz factor and isotropic radiation energy of X-ray flares (rest frame 0.3-10 keV, 
triangles) and prompt GRBs (solid dots) for the ISM case. The best fit power-law index of the 
correlation for prompt GRBs is 0.25. 




















Fig. 5.— Lorentz factor and isotropic radiation energy of X-ray flares (rest frame 0.3-10 keV, 
triangles) and prompt GRBs (solid dots) for the wind case. The best fit power-law index of the 
correlation for prompt GRBs is 0.33. 



Fig. 6.— Prompt (solid dots) and X-ray flare (triangles) Lorentz factors of 5 GRBs. 









































Table 1. Properties of GRB X-ray Flares 


GRB 

z 

T r i S e 

‘decay 

S x 

t 90,x 

E x ,iso, 50 

qISM 

gWind 

Lower rf M 

Lower rf ,w 

Upper F x 



( S) 

(s) 

(10 -7 erg cm -2 ) 

(s) 

n0 50 erg) 

(rad) 

(rad) 





050416 A 

0.650 

2 . 2 E 5 ± 2 . 8 E 6 

2 . 8 E 5 ± 5 . 8 E 5 

0.34 ± 0.53 

5 . 0 E 5 

0.4 ± 0.6 

0.026 ± 0.005 

0.071 ± 0.009 

43.7 ± 281.9 

15.9 ± 102.7 

14.6 ± 5.7 

050802 

1.710 

123.0 ± 0 

22.0 ± 0 

0.20 ± 0.30 

145.0 

1.5 ± 2.3 

0.028 ± 0.002 

0.041 ± IE -3 

14.9 ± 0 

10.3 ± 0 

157.1 ± 58.9 

050814 ( 1 ) 

5.300 

217.0 ± 0 

624.0 ± 0 

0.04 ± 0 

841.0 

2.0 ±0 

0.046 ± 0.003 

0.044 ± 0.002 

36.6 ± 0 

38.7 ± 0 

108.5 ± 0 

050814 ( 2 ) 

5.300 

505.0 ± 0 

439.0 ± 0 

0.05 ± 0 

944.0 

2.5 ±0 

0.046 ± 0.003 

0.044 ± 0.002 

20.1 ±0 

21.3 ± 0 

111.4 ±0 

050820 A 

2.620 

34.0 ± 0 

148.0 ± 0 

6.89 ± 0 

182.0 

110.5 ±0 

0.169 ± 0.007 

0.106 ± 0.003 

12.3 ± 0 

19.7 ± 0 

434.5 ± 0 

050904 ( 1 ) 

6.290 

120.0 ± 0 

107.0 ± 0 

2.51 ± 0 

227.0 

160.3 ± 0 

0.045 ± 0.007 

0.028 ± 0.003 

20.8 ± 0 

33.2 ± 0 

451.3 ± 0 

050904 ( 2 ) 

6.290 

96.0 ± 0 

188.0 ± 0 

0.27 ± 0 

284.0 

17.2 ± 0 

0.045 ± 0.007 

0.028 ± 0.003 

30.9 ± 0 

49.2 ± 0 

244.4 ± 0 

050904 ( 3 ) 

6.290 

86.0 ± 0 

108.0 ± 0 

0.11 ±0 

194.0 

7.0 ±0 

0.045 ± 0.007 

0.028 ± 0.003 

24.7 ± 0 

39.4 ± 0 

214.7 ± 0 

050904 ( 4 ) 

6.290 

1680.0 ± 0 

2236.0 ± 0 

0.88 ± 0 

3916.0 

56.2 ± 0 

0.045 ± 0.007 

0.028 ± 0.003 

25.4 ± 0 

40.6 ± 0 

170.4 ± 0 

050904 ( 5 ) 

6.290 

1176.0 ± 0 

7537.0 ± 0 

0.95 ± 0 

8713.0 

60.7 ± 0 

0.045 ± 0.007 

0.028 ± 0.003 

55.9 ± 0 

89.1 ± 0 

142.2 ± 0 

050904 ( 6 ) 

6.290 

5773.0 ± 0 

14457.0 ± 0 

0.60 ± 0 

20230.0 

38.3 ± 0 

0.045 ± 0.007 

0.028 ± 0.003 

34.9 ± 0 

55.7 ± 0 

102.7 ± 0 

050904 ( 7 ) 

6.290 

3774.0 ± 0 

1586.0 ± 0 

0.40 ± 0 

5360.0 

25.6 ± 0 

0.045 ± 0.007 

0.028 ± 0.003 

14.3 ± 0 

22.8 ± 0 

129.3 ± 0 

051016 B 

0.940 

109.0 ± 0 

1457.0 ± 0 

0.18 ± 0 

1566.0 

0.4 ±0 

0.162 ± 0.039 

0.242 ± 0.039 

22.6 ± 0 

15.1 ± 0 

63.3 ± 0 

060115 

3.530 

19.1 ± 10.2 

64.3 ± 9.3 

0.17 ± 0.07 

83.4 

4.3 ± 1.8 

0.044 ± 0.007 

0.047 ± 0.005 

41.8 ± 13.4 

38.7 ± 11.5 

234.5 ± 24.1 

060124 ( 1 ) 

2.300 

291.0 ±0 

70.0 ± 0 

27.13 ± 0 

361.0 

347.5 ± 0 

0.071 ± 0.005 

0.092 ± 0.005 

6.9 ±0 

5.3 ±0 

487.6 ± 0 

060124 ( 2 ) 

2.300 

50.0 ± 0 

313.0 ±0 

12.40 ± 0 

363.0 

158.8 ±0 

0.071 ± 0.005 

0.092 ± 0.005 

35.2 ± 0 

27.2 ± 0 

400.3 ± 0 

060210 ( 1 ) 

3.910 

23.1 ± 1.5 

37.6 ± 2.1 

2.30 ± 0.40 

60.7 

68.4 ± 11.9 

0.028 ± 0.003 

0.026 ± 0.002 

45.9 ± 6.1 

49.4 ± 4.6 

507.1 ± 22.0 

060210 ( 2 ) 

3.910 

14.7 ± 1.8 

45.5 ± 3.2 

1.20 ± 0.24 

60.2 

35.7 ± 7.1 

0.028 ± 0.003 

0.026 ± 0.002 

63.3 ± 9.1 

68.1 ± 7.4 

431.9 ± 21.6 

060418 

1.490 

5.6 ± 0.5 

19.5 ± 0.5 

4.80 ± 0.40 

25.0 

27.0 ± 2.3 

0.025 ± 0.005 

0.030 ± 0.004 

75.6 ± 16.6 

62.2 ± 9.3 

501.9 ± 10.5 

060526 ( 1 ) 

3.220 

11.8 ± 0.6 

12.7 ± 1.6 

3.20 ± 3.90 

24.5 

69.6 ± 84.9 

0.0850 ± 7 . 891 E -4 

0.0780 ± 4 . 86 E -4 

12.3 ± 0.8 

13.3 ± 0.9 

639.1 ± 194.7 

060526 ( 2 ) 

3.220 

9.8 ± 1.6 

27.9 ± 8.6 

2.50 ± 1.10 

37.8 

54.4 ± 23.9 

0.0850 ± 7 . 891 E -4 

0.0780 ± 4 . 86 E -4 

20.0 ± 3.5 

21.6 ± 3.8 

539.1 ± 59.3 

060526 ( 3 ) 

3.220 

15.2 ± 1.2 

34.6 ± 4.3 

2.70 ± 0.63 

49.9 

58.7 ± 13.7 

0.0850 ± 7 . 891 E -4 

0.0780 ± 4 . 86 E -4 

17.9 ± 1.3 

19.3 ± 1.4 

512.7 ± 29.9 

060526 ( 4 ) 

3.220 

10.2 ± 4.2 

61.7 ± 3.7 

1.30 ± 0.34 

71.9 

28.3 ± 7.4 

0.0850 ± 7 . 891 E -4 

0.0780 ± 4 . 86 E -4 

29.1 ± 6.1 

31.5 ± 6.6 

389.8 ± 25.5 

060707 

3.425 

8.1 ± 8.0 

26.5 ± 10.0 

0.07 ± 0.06 

34.6 

1.7 ± 1.4 

0.149 ± 0.024 

0.110 ± 0.012 

12.2 ± 6.7 

16.5 ± 8.9 

231.3 ± 49.6 

060714 ( 1 ) 

2.711 

8.1 ± 3.0 

41.6 ± 2.0 

3.30 ± 0.43 

49.6 

53.9 ± 7.0 

0.0270 ± 8 . 414 E -4 

0.0330 ± 6 . 938 E -4 

84.2 ± 15.9 

68.0 ± 12.8 

502.7 ± 16.4 

060714 ( 2 ) 

2.711 

4.0 ± 2.1 

4.2 ± 2.1 

0.38 ± 0.13 

8.2 

6.2 ± 2.1 

0.0270 ± 8 . 414 E -4 

0.0330 ± 6 . 938 E -4 

38.1 ± 13.9 

30.8 ± 11.2 

459.2 ± 39.3 

060714 ( 3 ) 

2.711 

4.4 ± 1.1 

9.6 ± 1.0 

0.95 ± 0.24 

14.0 

15.5 ± 3.9 

0.0270 ± 8 . 414 E -4 

0.0330 ± 6 . 938 E -4 

54.9 ± 7.6 

44.4 ± 6.1 

505.2 ± 31.9 

060714 ( 4 ) 

2.711 

8.3 ± 0.6 

12.8 ± 0.5 

1.10 ± 0.20 

21.1 

18.0 ± 3.3 

0.0270 ± 8 . 414 E -4 

0.0330 ± 6 . 938 E -4 

46.1 ± 2.4 

37.3 ± 1.7 

473.0 ± 21.5 

060729 

0.540 

9.6 ± 1.0 

34.1 ± 2.0 

9.30 ± 1.00 

43.7 

7.0 ± 0.8 

0.418 ± 0.037 

0.350 ± 0.020 

4.5 ± 0.5 

5.4 ± 0.5 

311.3 ± 8.4 

060814 

0.840 

9.3 ± 1.3 

29.7 ± 2.3 

3.10 ± 0.48 

39.0 

5.8 ± 0.9 

0.064 ± 0.006 

0.061 ± 0.004 

28.1 ± 3.5 

29.2 ± 3.0 

305.2 ± 11.8 

060906 ( 1 ) 

3.690 

349.0 ± 397.0 

1110.0 ± 1420.0 

0.05 ± 0.08 

1460.0 

1.4 ± 2.3 

0.027 ± 0.002 

0.032 ± IE -3 

66.9 ± 57.4 

56.0 ± 48.0 

86.8 ± 34.7 

060906 ( 2 ) 

3.690 

774.0 ± 987.0 

5030.0 ± 1640.0 

0.12 ± 0.06 

5810.0 

3.4 ± 1.7 

0.027 ± 0.002 

0.032 ± IE -3 

95.6 ± 63.3 

80.1 ± 52.9 

76.5 ± 9.6 

070306 

1.500 

8.7 ± 1.4 

34.7 ± 3.3 

2.10 ± 0.35 

43.4 

12.0 ± 2.0 

0.079 ± 0.019 

0.071 ± 0.011 

25.3 ± 6.6 

28.3 ± 5.3 

356.8 ± 14.9 

070318 ( 1 ) 

0.840 

9.5 ± 16.2 

19.4 ± 26.1 

0.10 ± 0.11 

28.9 

0.2 ± 0.2 

0.163 ± 0.011 

0.157 ± 0.007 

8.8 ± 9.6 

9.1 ± 9.9 

139.4 ± 38.3 

070318 ( 2 ) 

0.840 

44.9 ± 8.3 

103.0 ± 12.6 

1.30 ± 0.19 

147.9 

2.4 ± 0.4 

0.163 ± 0.011 

0.157 ± 0.007 

9.3 ± 1.2 

9.6 ± 1.1 

176.0 ± 6.4 

07072 IB ( 1 ) 

3.626 

0.9 ± 145.7 

120.8 ± 149.1 

0.60 ± 0.87 

121.7 

15.8 ± 23.0 

0.0210 ± 7 . 033 E -4 

0.0230 ± 5 . 182 E -4 

561.6 ± 45459.9 

508.2 ± 41133.9 

295.7 ± 107.2 

07072 IB ( 2 ) 

3.626 

4.4 ± 4.4 

7.3 ± 5.5 

0.26 ± 0.18 

11.7 

6.9 ± 4.8 

0.0210 ± 7 . 033 E -4 

0.0230 ± 5 . 182 E -4 

62.4 ± 39.1 

56.5 ± 35.4 

430.8 ± 74.6 

07072 IB ( 3 ) 

3.626 

8.9 ± 4.0 

18.3 ± 6.3 

0.19 ± 0.07 

27.2 

5.0 ± 1.8 

0.0210 ± 7 . 033 E -4 

0.0230 ± 5 . 182 E -4 

69.5 ± 19.8 

62.9 ± 17.9 

322.6 ± 29.7 

07072 IB ( 4 ) 

3.626 

62.2 ± 64.3 

179.8 ± 151.2 

0.22 ± 0.18 

242.0 

5.8 ± 4.8 

0.0210 ± 7 . 033 E -4 

0.0230 ± 5 . 182 E -4 

82.4 ± 55.0 

74.6 ± 49.7 

193.7 ± 39.6 

071031 ( 1 ) 

2.690 

25.0 ± 3.5 

36.2 ± 2.6 

2.80 ± 1.10 

61.2 

45.2 ± 17.8 

0.061 ± 0.011 

0.069 ± 0.009 

19.8 ± 4.0 

17.3 ± 2.5 

456.4 ± 44.8 

071031 ( 2 ) 

2.690 

9.3 ± 1.9 

34.7 ± 5.0 

0.80 ± 0.19 

44.0 

12.9 ± 3.1 

0.061 ± 0.011 

0.069 ± 0.009 

31.8 ± 7.1 

27.8 ± 4.9 

362.4 ± 21.5 

071031 ( 3 ) 

2.690 

9.8 ± 2.4 

42.1 ± 6.3 

0.54 ± 0.13 

51.9 

8.7 ± 2.1 

0.061 ± 0.011 

0.069 ± 0.009 

34.1 ± 8.0 

29.8 ± 5.7 

315.2 ± 19.0 

071031 ( 4 ) 

2.690 

65.5 ± 4.4 

210.6 ± 9.0 

1.90 ± 0.19 

276.1 

30.7 ± 3.1 

0.061 ± 0.011 

0.069 ± 0.009 

29.5 ± 5.6 

25.8 ± 3.3 

284.3 ± 7.1 






